Sommerfeld Fine-Structure Formula for Two-Body Atoms 
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For relativistic atomic two-body systems such as the hydrogen atom, positronium, and muon- 
proton bound states, a two-body generalisation ol the single-particle Sommerfeld fine-structure 
formula for the relativistic bound-state energies is found. The two-body Sommerfeld bound-state 
energy formula is obtained from a two-body wave equation which is physically correct to order 
(Za) 4 . The two-body Sommerfeld formula makes two predictions in order (Za) 6 for every bound 
state and every mass ratio. With N the Bohr quantum number: (a) The coefficient of the (Za) 6 /N 6 
energy term has a specified value which depends only on the masses of the bound particles, not on 
angular quantum numbers; (b) The coefficient of the (Za) 6 /N 4 energy term is a specified multiple 
of the square of the coefficient of the (Za) 4 /N 3 energy term. Both these predictions are verified in 
positronium by previous calculations to order (Za) 6 which used second-order perturbation theory. 
They are also correct in the Coulomb-Dirac limit. The effect of the two-body Sommerfeld formula 
on calculations of muon-proton bound- state energies is examined. 
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I. INTRODUCTION 

The lack of an analytically solvable relativistic wave 
equation for two-body atomic systems has compelled 
physicists to use second-order perturbation theory in cal- 
culating energy levels to order (Za) e in systems such 
as positronium 0, Q . Several two-body wave equations 
or prescriptions (reviewed in ref. |3|), starting with 
the Breit equation Q in 1929, have correctly described 
the energy levels of two-body atomic systems to order 
(Za) 4 using first-order perturbation theory. But none 
has been solvable analytically, so none has been able to 
be a starting-point for the calculation of atomic energy 
levels at order (Za) 6 without the use of second-order per- 
turbation theory. 

For the rest of this paper the symbol a will be a syn- 
onym for Za. No radiative corrections are considered. 

To explain what is found in this paper and what re- 
mains to be found, we remind the reader that the ana- 
lytic solution [B| to the Coulomb Dirac equation in 1928 
gave the one-particle Sommerfeld fine-structure formula 
for the bound-state energies: 
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Here m is the mass of the electron, a is the fine-structure 
constant, n is the radial quantum number and enirac is 
determined by the eigenvalue of a finite-dimensional an- 
gular equation 



(eoirac + 1) 2 £ = [(1 + <r ■ Lf - a 2 ]£. 



(2) 



Because equation ([2]) for the angular parameter e Dirac 
is exactly solvable, the bound-state energies predicted 
by equation ([!} can be evaluated by simple algebra not 



only to order a 4 , but to order a e and even a 8 . The 
Coulomb-Dirac problem does not need second- and third- 
order perturbation theory. 

If a two-body wave equation accurate to order (Za) 4 
were to be solvable analytically like the one-particle 
Coulomb Dirac equation, no second-order perturbation 
theory would be needed to obtain two-body energy levels 
to order (Za) 6 . The other terms needed to complete the 
physical calculation to order (Za) 6 would only need to 
be evaluated in first-order perturbation theory. 



II. MAIN RESULT 

In this paper we will find a two-particle counterpart to 
the Sommerfeld fine-structure formula {1}. It is 
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Here M is the mass of the other particle. It is easy to see 
that E — M reduces to the one-particle Sommerfeld fine- 
structure formula (QJ when M — > oo. The formula Q 
for the bound-state energies is obtained from a two-body 
relativistic atomic wave equation which was derived in 
ref. [3| and is repeated below. 

In this work we will not be able to derive an angu- 
lar equation for the angular parameter e which is a two- 
particle counterpart to equation Therefore we have 
only half the analytic solution. 

Nevertheless, from the two-particle Sommerfeld for- 
mula ([3]) alone, assuming that e exists without knowing 
it, we have found two predictions for order a 6 bound- 
state energy terms for any masses, which are verified by 
previous calculations for positronium which had to use 
second-order perturbation theory. These results raise the 
hope that it may eventually be possible to find a full 
analytic solution to a relativistic atomic two-body wave 
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equation which would eliminate the need for second-order 
perturbation calculations to obtain energies to order a 6 . 



III. TWO-BODY WAVE EQUATION 

The local relativistic atomic two-body wave equa- 
tion © below was derived in ref. Q from the Bethe- 
Salpeter equation for two spin-1/2 constituent point 
particles bound by a singlc-photon-exchange kernel in 
the Coulomb gauge. The derivation used a simple 
quasi-potential approximation [(| @] with its associated 
Blankenbecler-Sugar correction series Q . 

The masses of the bound particles are denoted by 
m, M . The bound-state energy E is parametrised by a 
quantity /3 as follows: 



E = yjm 2 ~/3 2 + ^AP - P 2 



(4) 



The particles' individual bound-state energies also occur: 



t = yjm 2 - /3 2 , T = y/M 2 - /? 2 



(5) 



The Pauli matrices <x and the Dirac matrices 7, 7 refer 
to the particle of mass m, while £, T and T° refer to the 
particle of mass M. The operator p is — zV, where V 
refers to the relative position coordinate r. Also r = |r| 
and f = r/r. 

The relativistic bound-state wave equation in the 
ccntrc-of-mass system is 

[p 2 + /3 2 Mr) 
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The eigenvalue is /3 2 , not _E. The constant (3 2 is substi- 
tuted into the square roots in equation (j4|) to obtain the 
energy E. There are no non-local operators of the form 
\/p 2 + m 2 , so no kinetic-energy correction terms p 4 , p 6 
appear in perturbation theory. 

The first and second terms in the curly brackets are the 
standard binding potential and Breit interaction. The 
third term is an effective term equal to 1/2E times the 
square of the binding potential. This is a consequence 
of the Blankenbecler-Sugar correction formalism, which 
ensures that the relativistic energies are correct The 
final term in the curly brackets, which was not included 
in ref. [3j, contains an anomalous magnetic moment for 
the particle of mass M. It is derived by adding a term 
i(g — 2)(k x S)/4M to the Dirac matrix T in the ver- 
tex function of the particle of mass M in the originating 
Bcthe-Salpctcr equation. The term is included here so 
that the particle of mass M may represent a proton. 

As discussed in ref. @, the wave equation © 
is expected to be valid for r > a/2E, because the 



Blankenbecler-Sugar correction formalism results in a se- 
ries in powers of (— a/r)/2E. In addition, for the hydro- 
gen atom and muonic hydrogen [13, EH, a/2E is about 
0.001 fm. The mean-square charge radius of the proton 
is about 1 fm. Therefore except perhaps for positron- 
ium the validity of the two-body wave equation ([6]) is 
expected to end well before r descends to a/2E. 

In the limit M 00, dividing through by the first 
right-hand operator in square brackets, it is easy to see 
that equation ([6]) reduces to the one-body Coulomb Dirac 
equation 



[m + 7 • p — 7°i] ^(r) 
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which gives the energies ([T]) and the angular equation ([2j . 

To verify the correctness of the wave equation © we 
give its bound-state energies to order a 4 . With J (often 
written F) denoting the total spin of the bound state, 
notation used is: 
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From first-order perturbation theory with iV the Bohr 
quantum number the bound-state energies of the wave 
equation ([6]) are 



E = m + M - 
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in which 
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L = J + 1 
L = J — 1 
L = J, S « 1 
L = J, S w 

(11) 



$ = 



1-2x3/3 L = 0,5 = 1 



1 + 2xg 



L = 0,5 = 



(12) 



The symbols 5 ~ 1, 5~0 for L = J stand for the state 
in which S is predominantly 1 or 0, respectively. 

Examining the results above, we find that in the limit 
M — > 00 the energies of the Dirac Coulomb equation 
are recovered. For m -C M it is easy to verify that the 
hyperfine splittings are correct to order m/M ([Hj], Sec. 
22). With m — M and g = 2 we find the standard 
positronium energies without the annihilation term ( [l2j , 
Sec. 23). These results verify that to order a 4 , the two- 
body wave equation ^ gives physically correct results. 
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IV. DERIVATION OF THE TWO-BODY 
ENERGY FORMULA 

Following conventional treatments of the Coulomb 
Schrodinger equation and the Coulomb Dirac equation, 
we substitute 



(13) 



that the wavefunction remain Unite at large r, we cannot 
use the behaviour of eqn. © as r — >• to determine e. 

However, the assumption that the quantity e does exist 
will lead below to some new predictions which arc verified 
in known physical systems such as the hydrogen atom 
and positronium. The predictions apply to muon-proton 
bound states as well. 



into the wave equation 



Here the coefficients cu are 



16-dimensional vectors. The expansion is expected to be 
valid for r > a/225, or for r greater than the proton 
radius if a proton is present. One obtains a four-term 
recurrence relation for the coefficients aj. 

When the dominant terms acting on the large compo- 
nent of the wave function in the wave equation ([5]) are ex- 
amined, we see that they are the same as in the Coulomb 
Schrodinger equation. That means that if the series (TIB")) 
does not terminate for some j — n, the wavefunction will 
diverge as e +l3r for large r. So that the wavefunction (fT5j) 
will go to zero at infinity, in this paper we will assume 
that the series terminates at j = n. Then it is easy to 
show that with aj = for j > n, the recurrence relation 
implies this equation for a n : 



2(3(e + 1 + n)a v 



= fhM 
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containing the abbreviations 



m 



■ i/3-y ■ r + j°t, M = M + if3T-r + T°T. (15) 



V. EXPANSION OF E IN POWERS OF a 2 

In the bound-state energy formula ([3]) the quantity n + 
e+1 occurs. While in the Coulomb Schrodinger equation 
e is always L, in the Coulomb Dirac case equation @ 
sometimes gives eDirae ~ L, but sometimes eDirac ~ L— 1. 
In the latter case the ground state is known to be given 
by n = 1, not n — 0. Thus in either case, n + e Dirac + 1 is 
still approximately the Bohr quantum number N. In fact 
it is easy to find from @ that, with k = |L + cr/2\ + 1/2, 



n + EDirac + 1 = N - 



a 



2k 8k 3 16k 5 



(17) 



It is likely that the same thing happens for the two-body 
wave equation (j6|). Therefore we expand 



n + e+l = N + a 2 e 2 + a 4 e 4 + a 6 e 6 



(18) 



When the expansion (fT8|) is substituted into the two- 
body Sommerfeld formula ([3|) , the energy levels E can be 
expanded in powers of a 2 . Recalling that x = [i 2 /mM, 
one finds 



Equation (| 14[) shows that a n — rhMb. Substituting 
that back into (|14p puts a projection operator fhM on 
each side of the Coulomb and Breit terms; multiplying 
them out gives a scalar multiple of fhM again. One thus 
finds 



2/3(e + r. 
which is to say 
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Using equations @ and (O, equation (Tip]) gives the 
two-body Sommerfeld "fine-structure" bound-state en- 
ergy formula ([3]). 

The parameter e, which depends on the angular quan- 
tum numbers, cannot be found from the two-body wave 
equation (J6j) as easily as e Dirac is found from the Coulomb 
Dirac equation (|7|). The reason is that eoirac is deter- 
mined by the regular singularity of ([7]) as r — » 0. But 
for finite m and M, the wave equation ^ is expected to 
be valid only for r exceeding the small but finite value 
a/2E, or exceeding the radius of a proton if a proton is 
one of the bound particles. While the two-body Sommer- 
feld energy formula (j3]) has come from the requirement 



E = m+M+a 2 nC 2 +CY 4 nC4+a 6 fiC 6 +a 8 nC s 



in which 



C 2 = - 

c 4 = 
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(20) 
(21) 



5 - 3z + x 2 (3-z)e 2 3e| e 4 , 00 , 

°6 TTHTFk 771^ 77777 + T77 \2l) 
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16iV 6 2N 5 2N 4 N 3 

35 - 29a; + 18a; 2 - 5a; 3 



128A^ 8 



3 [5 - 3a; + a; 2 ] e 2 5(3 - x)e 



8N 7 



(23) 



4e 2 - (3 - x)e 4 3e 2 e 4 e 6 



2N 5 



N 4 N 3 



In the Coulomb-Dirac limit, fi = m and x = 0. In this 
limit, using the expansion of eoirac from eqn. (|17p . eqns. 
P^)l to ([2^|l reproduce the known expansion of i?Dirac to 
order a 8 . For general masses m,M we now examine the 
expansion (TT9)) and the coefficients ([20)) to (|23l) to see 
whether first, they agree with the known results to order 
a , and second, to see whether they predict anything not 
known yet. 
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VI. a 2 AND q 4 TERMS 

The coefficient Ci gives the Bohr levels correctly. 

Next, referring to the energy levels (fTU)l given above, 
even the coefficient C4 suggests the correctness of the 
two-body Sommcrfeld formula ©. Eqn. (|2"T|) agrees 
that there are a 4 /N 4 and a 4 /N 3 terms only. Moreover, 
as eqn. (fT0|) shows, C4 gives the a 4 /N 4 term correctly 
for all masses, verifying that it is independent of the un- 
known value of e. Finally, from the first-order perturba- 
tion theory result (TIU)) we can read off the value of £2: 

ea = ~[(l-*zo)6 + *io*] (24) 
with 9 and $ given by §TQ and (Tig) . 

VII. PREDICTIONS OF a 6 TERMS 

The a 6 terms in (l2"2"j) will tell whether the two-body 
Sommcrfeld energy formula ([3]) contains any information 
previously obtainable only from using second-order per- 
turbation theory, and whether it predicts anything as yet 
uncalculated. Note that the predictions described hold 
automatically in the Coulomb-Dirac limit M — > 00. 

A. a 6 /N e > Terms 

The coefficient of the a 6 /j,/N 6 term given by (|2"21 does 
not need knowledge of the angular eigenvalue e. It is 
predicted to have the same value — [5 — 3x + £ 2 ]/16 for 
every state, depending only on the masses. Therefore 
for positronium, for which [i = m/2 and x — 1/4, eqn. 
(|2"2")l predicts that the coefficient of a 6 m/N 6 is always 
—69/512. In fact, precisely this value has been found for 
all P-states for all S-states 0], and in Zatorski's re- 
cent calculations for all L > 2 states [l|. To contrast the 
methods, we note that in ref. [J the coefficient -69/512 
is the sum of five terms, three from second-order pertur- 
bation theory (eqns. (159), (160) and (162) of ref. [l[) 
and two first-order terms, one of which contains the ex- 
pectation value of the p 6 /to 5 kinetic-energy correction 
which is absent in our formalism (ref. [l[ eqns. (91) and 
(122)). 

In addition, Pachucki |14j has found the same function 
— [5 — 3x + x 2 ]/16 for the coefficient of the a 6 n/N 6 energy 
term as in (|22[) for arbitrary masses — but only for S- 
states. Accordingly this coefficient remains to be tested 
for L > 1 when to ^ M. 

B. a 6 /N 4 Terms 

For L > 1, in ref. [l[ the a 6 to/A 4 terms of the positro- 
nium energy levels come from the sum of six second-order 
perturbation terms (eqns. (153), (158), (164), (171), 



(174) and (177)) and no first-order terms. Since the two- 
body Sommerfeld formula ^ is a candidate to replace 
the results of second-order perturbation theory once e 
is fully determined, these calculated a 6 /N 4 terms are a 
possible second test of the Sommerfeld formula ([3]). 

In positronium, from (|2ip the a 4 /N 3 energy term is 
(e 2 /2)a 4 m/iV 3 , and from the a 6 /N 4 energy term 
is (— 3e2/4)a 6 m/iV 4 . Therefore it is predicted that the 
coefficient of the a 6 m/N 4 term of positronium should be 
—3 times the square of the coefficient of the a 4 m/N 3 
term. 

As an example for L — J + 1, with m = M and g = 2, 
from (fTU)) and (fTTj) above the coefficient of a 4 m/N 3 in 
positronium is 

-(4L 2 + L- 1) 
4L(2i- 1)(2L+ 1)' 

In eqn. (211) of ref. |l| the calculated coefficient of the 
a 6 m/N 4 term for L = J + 1 is given from the sum of the 
six second-order terms as 

3 - 6L - 21L 2 + 24L 3 + 48L 4 
16L 2 (2L- 1) 2 (2L+ l) 2 

It can be seen that the latter expression is indeed —3 
times the square of the former. This verifies the predic- 
tion of the two-body Sommcrfeld formula for positron- 
ium L = J+ 1 states. The prediction also gives the other 
a 6 m/N 4 terms of positronium energies correctly (eqns. 
(207), (215) and (219) of ref. [3). 

For the 1 5o and 3 Si states of positronium the coeffi- 
cients of the a 6 m/N term contain other contributions 
(see ref. [2[) and no prediction can be made. For L > 1 
and to 7^ M the prediction for the a 6 /N 4 term remains 
to be tested. 



C. q 6 /A 5 and a 6 /N 3 Terms 

While (|2"2"1) predicts the coefficient of a 6 ^/A 5 to be 
— (3 — x)e2/2, in which £2 is already known from (|24[) . 
we cannot tell whether this coefficient is confirmed by 
previous positronium calculations because many terms 
of order a 6 contain 1/N 5 [HQ. 

The final coefficent £4 of a 6 fi/N 3 in (|22p cannot be 
evaluated and tested until an angular equation for e anal- 
ogous to (T2J) is found from which £4 can be calculated. 
Note, however, that since L = J ± 1 mixing is prevented 
in order a 4 because (L = J — l\l /r 3 \L = J +1) =0, this 
mixing does not occur in £2 either, because of eqn. (|24[) . 
Consequently, such mixing can only occur in the £4 term. 
Therefore eqn. (f2"2")) predicts that in order a 6 , L = J± 1 
mixing should occur in the a 6 /N 3 term only. 
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VIII. OTHER REMARKS 

A. m < M 

For m <C M and e = eDirac the two-body Sommerfeld 
formula reproduces a result in the review of hydrogen- 
like atoms by Eides, Grotch and Shelyuto [15} . The two- 
body Sommerfeld formula ([3} can be written, once again 
recalling that x = fj,/(m + M), as 



E = (m + M)VT+2xh 



with 



1 



(25) 



(26) 



(n+e+l) 2 



For m <C M, x is small. Equation (1231) expanded in 
powers of x up to order x 2 is 



E = (m + M)+nh 



1 /x 



2m + M 



(27) 



With e = eDirac this is the same expression that appears 
in eqn. (38) of ref. [l5j]. Expansion of expression ([27} 
to order a 6 reproduces the expression ([22} for Cg except 
for the x 2 term, which is small for the hydrogen atom. 
This agreement represents another confirmation of the 
two-body Sommerfeld formula ([3]). 



Then the Blankenbecler-Sugar correction formalism Q 
can be used to evaluate the remaining corrections. 

Since for real two-body atoms the wave equation © 
has not revealed a way to calculate e yet, we report on 
a complete implementation of part (i) for a far simpler 
model. 

Ref. [l6| used a Bethe-Salpeter equation for scalar 
constituent particles of masses m, M bound by a scalar- 
Coulomb kernel —AraMa/r. Following through the for- 
malism of ref. [3} gave a two-body wave equation with a 
solution which could be evaluated fully. The wave equa- 
tion had a regular singularity as r — > which gave an 
algebraic equation for e not unlike eqn. ([2} of the present 
paper. The solution series terminated at j = n for ap- 
propriate energies. The resultant two-body Sommerfeld 
formula was quite similar to eqn. Q of the present pa- 
per. It gave energies to order a 6 since e was known. The 
Blankenbecler-Sugar correction formula was used to find 
the corrections to the approximations which led to the 
wave equation. Evaluation of their expectation values 
gave a 6 terms for all L, as well as a 5 and a 6 log a terms 
for L = 0. Thus the calculation of the bound-state en- 
ergies of the model scalar Bethe-Salpeter equation was 
completed to order (Za) e without using second-order 
perturbation theory. 

It is possible that a similar implementation for real 
two-body atoms using a Coulomb-gauge kernel would re- 
quire a return to the derivation in ref. Q, perhaps ob- 
taining a wave equation sightly different to © near the 
origin. 



B. 



(Za) 



The expansion of the bound-state energy given by 
equations (fT9|) to ([23} is carried out to order a 8 , meaning 
(Za) 8 . As a single example, in the unlikely event that 
any energy levels of positronium are calculated to order 
(Za) 8 for some angular quantum numbers and all Bohr 
quantum numbers N, eqn. (|23[) predicts that the coeffi- 
cient of (Za) s m/N s for every state of positronium will 
be 1843/16384. 



D. Muonic Hydrogen 

For muon-proton bound states we have compared the 
energies calculated by first-order perturbation theory 
((TO} to an evaluation using the exact two-body Sommer- 
feld formula ([3} in which, given the value of n + e + 1, 
recoil effects appear to be exact. 

Substituting the approximation 



C. Evaluation from Bethe-Salpeter Kernel 

The wave equation ([6} was derived from a Bethe- 
Salpeter equation whose kernel consists of the exchange 
of a single Coulomb-gauge photon. To evaluate physical 
energies to order (Za) 6 fully, the procedure would be (i) 
to evaluate the energies of this binding Bethe-Salpeter 
equation fully to order (Za) 6 , then (ii) to evaluate the 
energy contributions of other relevant kernels, including 
crossed diagrams, all of which would need only first-order 
expectation values to obtain energies to order (Za) 6 . 

In part (i) the first task is evidently to find and solve a 
two-body angular equation similar to ([5} for the unknown 
angular parameter e, so that the two-body Sommerfeld 
formula @ can be evaluated to any desired accuracy. 



n + e + l = N + a 2 e 2 

into the two-body Sommerfeld formula @ with e 2 given 
by the first-order perturbation value (|2~4"} . it was found 
that the energy of the nominal lSfy^ state evaluated by 
the two-body Sommerfeld formula was 0.005 meV lower 
than when evaluated by first-order perturbation theory. 
The other IS state had a smaller difference. The energy 
of the nominal 2S^° state was 0.0008 meV lower when 
calculated by the Sommerfeld formula than by perturba- 
tion theory. The other N = 2 states had a smaller differ- 
ence. Thus for muonic hydrogen the recoil correction to 
first-order perturbation theory given by the Sommerfeld 
formula appears to be small for the 2P and 2S states of 
current interest [10, [ill • 
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IX. CONCLUSION 

Verification of the two (Za) 6 predictions, especially 
for the (Za) 6 /N 4 term which was quite unexpected, in- 
dicates two things, (i) The two-body Sommerfeld energy 
formula (|3|) seems to be true, (ii) The angular eigenvalue 
e apparently exists. We do not know yet whether the 
two-body atomic wave equation ([6|), valid as written for 
r > a/2E, determines a hnite-dimensional angular equa- 
tion for e analogous to the one-particle angular equation 
© for 

It is expected that any complete analytic solution to 
the relativistic two-body atomic bound-state problem 
will give a two-body Sommerfeld formula for the bound- 
state energies which would give the (Za) 4 levels com- 
pletely, and which would also give the contribution of 
order (Za) 6 which previously has had to be calculated 
from second-order perturbation theory on the Breit in- 



teraction. To obtain complete physical results to order 
(Za) 6 , as was done for the simple scalar Bethe-Salpeter 
equation in ref. [16| , it is necessary first to go back to the 
originating Bethe-Salpeter equation and evaluate its en- 
ergies to order (Za) 6 . This is done by letting the formal- 
ism given in [3| determine what small additional terms 
are present beyond those contained in the Sommerfeld 
formula. The expectation values of these terms are eval- 
uated. 

More than 80 years have passed since the analytic solu- 
tion to the Coulomb Dirac equation gave the Sommerfeld 
fine-structure formula ([1]). The existence of the two-body 
Sommerfeld formula (J3)), even though its angular eigen- 
value e remains unknown for now, suggests that a wave 
equation for two-body atoms, perhaps somewhat similar 
to our equation (JSJ, may eventually be found and solved 
analytically. 



[1] J. Zatorski, Phys. Rev. A 78, 032103 (2008). 

[2] A. Czarnecki, K. Melnikov and A. Yelkhovski, Phys. Rev. 

A 59, 4316 (1999). 
[3] J.H. Connell, Phys. Rev. D 43, 1393 (1991). 
[4] G. Breit, Phys. Rev. 34, 553 (1929). 
[5] W. Gordon, Z. Physik 48, 11 (1928); C.G. Darwin, Proc. 

Roy. Soc. Ser. A 118, 654 (1928). 
[6] J. H. Connell, Thesis, University of Washington (Seattle, 

1967). 

[7] I.T. Todorov, Phys. Rev. D 3, 2351 (1971). 
[8] R. Blankenbecler and R. Sugar, Phys. Rev. 142, 1051 
(1966). 

[9] J.H. Connell, SLAC-PUB 5633 (1991). 
[10] A. Antognini, F. Kottmann, F. Biraben, P. Indelicato, F. 

Nez and R. Pohl, |arxiv . org/abs/1208 ■ 2637 1 
[11] R. Pohl, R. Gilman, G. A. Miller and K. Pachucki, Ann. 



Rev. Nucl. Part. Sci. (2013); [arxiv7o rg/abs/1301 . 0905 

[12] H.A. Bethe and E.E. Salpeter, Quantum Mechan- 
ics of One- and Two-Electron Atoms (Springer- 
Verlag/Academic, New York, 1957) [also published 
in Handbuch der Physik, Vol. XXXV (Springer, 
Berlin/Gottingen/Heidelberg, 1956)]. 

[13] LB. Khriplovich, A.I. Milstein and A.S. Yelkhovski, Zh. 
Exp. Teo. Fiz. 105, 299 (1994). This paper's results are 
quoted in refs. QJ and @]. 

[14] K. Pachucki, Phys. Rev. Letters 79, 4120 (1997). 

[15] M. I. Eides, H. Grotch & V. A. Shelyuto, Theory of 
Light Hydrogen-like Atoms (Springer Tracts in Mod- 
ern Physics 222, 2006); Phys. Rep. 342, 63 (2001); 
|arxiv . org/abs/hep-th/0002158"1 

[16] J. H. Connell, |arxiv7or g/abs/hep-th/0006082 



